ON THE LIFETIME OF QUASI-STATIONARY STATES IN 
NON-RELATIVISTIC QED 



DAVID HASLER, IRA HERBST, AND MATTHIAS HUBER 



Abstract. We consider resonances in the Pauli-Fierz model of non-relati- 
vistic QED. We use and slightly modify the analysis developed by Bach, 
Frohlich and Sigal [3j [4] to obtain an upper and lower bound on the lifetime 
of quasi-stationary states. 



1. Introduction and Main Result 

Spectral properties of models of non-relativistic QED were investigated by 
Bach, Frohlich, Sigal and Softer [2 O HI [5] and by many others. Bach, Frohlich, 
and Sigal [4] proved, among other things, an upper bound on the lifetime of quasi- 
stationary states. 

We show an upper and lower bound on the lifetime of quasi-stationary states. 
We heavily rely on the analysis developed in [3j [4] , but choose a different contour 
of integration and make use of an additional cancellation of terms. Moreover, we 
neither require a non-degeneracy assumption nor a spectral cutoff. However, we 
do not provide time dependent estimates on the remainder term and there are 
no photons in our quasi-stationary state. Estimates similar to ours were obtained 
before by different authors for other models, see e.g. [T6l[i7] . 

In order to be self-contained, we give all the necessary definitions for the 
model considered. For details, we refer the reader to [4]. We consider an atom in 
interaction with the second quantized electromagnetic field. The Hilbert space of 
the system is given by 

H := Hei ® T y 

where 

H e i :=A N L 2 l(R 3 x Z 2 )] N 
is the Hilbert space of N electrons with spin, and where 



T:= 05ArL 2 [(R 3 xZ 2 f 

jV=0 
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is the Fock space (with vacuum Q) of the quantized electromagnetic field, allowing 
two transverse polarizations of the photon. An and Sn are the projections onto 
the subspaces of functions anti-symmetric and symmetric, respectively, under a 
permutation of variables. Strictly speaking, we would have to take the physical 
units into account in the definition of these spaces. However, we refrain from 
doing so in order not to complicate the notation. The operator 
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l<i<j<A' 
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describes the electrons, and the operator for the total system is 



1 



N 

9 2m ^ 



(-ihV Xj - tA'Ax 3 ))] 2 ■■ +H' } 



N 

E 
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E ■ 

l<i<j<N 



where — eZ is the charge of the nucleus, e < the charge of the electron, c the 
velocity of light, h is Planck's constant, eo is the permittivity of the vacuum, m 
the mass of the electron, <jj is the Pauli matrix for the jth electron, and : • • • : 
denotes normal ordering. The kinetic energy of the photons is 



H> := he 



E 



dk\k\a>;{k)a^(k), 



where the a'*(k) and a'^k) are the usual creation and annihilation operators. The 
second quantized electromagnetic field is A' K , (x) := A! K , + A' K , where 



and 



A' K ,(x) + := 



A' A 



E 

M=l,2 



E 



dkK,'{\k\)\ 



dkn'{\k\) s 



2e c|fc|(27r) 



e'Jk) t 



3 f 



-ik-x /* 



(k). 



2e c|fc|(27r) 



3 s^k)e^a'^k) 



Here s'^(k), \x — 1,2, are the polarization vectors of the photon, depending only 
on the direction of k. Let us note that we use SI units here; for details about these 
operators, we refer the reader to [HI HO] . We set a := (Bohr radius), 

£ := and :— |2. Moreover, n'(r) := K(r^) is a cutoff function depending on 



the fine structure constant a 



47T€q^C 



k is a function, which is positive on [0, oo), 



satisfies k(t) — > 1 as r — > 0, and has an analytic continuation to a cone around the 
positive real axis which is bounded and decays faster than any inverse polynomial, 
e.g., k(t) :— e~ r . Following [3j, we scale the operator with the transformation 
Xj — > Qxj and k — > £ _1 fc. We denote the corresponding unitary transformation 
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by U. After this transformation the electron positions are measured in units of 
^ao, photon wave vectors in units of — , and energies in units of 4Ry, where the 

2 2 

Rydberg is Ry := a " c . The creation an annihilation operators transform as 

Ua'^k)!/- 1 = e /2 a„m, Ua';{k)U^ = e'^H^)- 
Moreover, we set 

£/.(*) :=<(r^), M=l,2. 

Accordingly, we obtain 

UH'gU- 1 = 2a 2 (mc 2 )H g , 
with H g := Hq + W g and Hq := H e i <8> 1/ + l e z <8> Hf, where 

N -Z 1 
H e i := -A 3N + 1—7 + ^ I— — |- 

3=1 l<i<i<JV R 



Here 



:= £ / d*|*|a*(*K(*) 



and the interaction is given by 

N 

W g : £{2a 3 ' (-*"V BJ ) + a 3 : A 2 (a^) : 

J'=l 

+ a 5/2 tTj • (V x A re )(aa:j)}, 

where the second quantized electromagnetic field is A K (ir) := A K (x) + + 
with 

4.0-0+ := Y I ^§ e( ,(*]e^(*). 

^=1,2^ V 4 ^ l fc l 

and 

As in [3], we set g := c?l 2 . Henceforth, we let the coupling constant g := 
a 3/2 > o be the perturbation parameter. We assume that the spectrum of H e i has 
the structure 

a(H e i) = {E ,E 1) ...}U[E,oc), 

where E := inf a ess (H e i) and Eq < E\ < . . . arc (at least two) eigenvalues (pos- 
sibly) accumulating at S. In the following, we will look at one (fixed) eigenvalue 
Ej of H e i with j > 1. For < e < 1/3 we set p a := g 2 ~ 2e , and A(5,e) := 
[Ej - 6/2, Ej + 5/2} + i[-g 2 ^, oo), where (5 := dist(Ej , a(H e i) \ {Ej}) > 0. We 
define the operators 

(1) H el {6) :=U el (6)H el U el {9)-\ H g (9) := U(9)H g U(9y 1 , 

W a (0) :=U{e)W g U{6)- 1 
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for real 6, where U (9) is the unitary group associated to the generator of dilations. 
It is defined in such a way that the space coordinates of the electrons are dilated 
as Xj i ► c Xj and the momentum coordinates of the photons as k i— > e k. It can 
be shown [4] Corollary 1.3, Corollary 1.4] that the operators defined in equation 
(fTJ) are analytic families for \9\ < 8q for some 9q > 0. We introduce the convention 
9 := I'd with d > 0. Moreover, U e i(9) is the above dilation acting on the electronic 
space only. 

We define (with r > small enough) P eUl {9) := -{2m)- 1 f, B ._ z ,(H e i {9) - 
z)~ 1 dz to be the projection onto the eigenspace corresponding to the eigenvalue 
Ei of H el (9) and set P e /,i(0) := 1 - P e i,i{9). Furthermore, we define P(9) := 
Pdj{6) (8 XH f < P0 and P(9) := 1 - P{9). We abbreviate P eM := P eM (0). 

Note that if we consider operators of the form PAP, where A is a closed 
operator and P a projection with DomA C RanP, then our notation does not 
distinguish between the operators PAP and PAP\^g, n p. It will be clear from the 
context, how the symbol PAP is to be understood. 

Following [1] , we make crucial use of the Feshbach operator 

(2) T P{e) (H g (9) - z) := P(9)(H g (9) - z)P{9) 

- P(9)W g (9)P(9)[P(9)(H g (9) - z)P(9T 1 P(9)W g (9)P(9). 

For the convenience of the reader, we summarize its most important properties 
including its existence in Appendix A. For details, we refer the reader to [31 Section 
IV] and [1] . It was shown in [31 0] that the Feshbach operator can be approximated 
in a sense to be shown using the operators 

(3) Zf(a):=]imJ2 [ dkP^w^k, //) 

X Pel,j[Pel,jH e l - Ej + \k\ - ie^PeljW^k, fl)P d J 

and 

(4) Zf(a):=J2 I ^PeLAl{k^)P el , 3 wfl{k^)P eUr 

^ h2 Jkm* 1*1 

Here the coupling functions WQ-[(fc,/i) and u4o(&,/i) will be needed later with 
9 ^ 0. Denoting the momentum of the jth electron by Pj, they are 

N 

(5) w { e }(k^) :=w®{k,») == £{2e- e Gf (fc, M ) ■ W + ^ • B«(M)}, 

j'=i 

where 

(6) Gf(fc, M ) := e ^ K X^ l) e lak - x e^(k) 
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and 

(7) flW(* lM ) := ae ~ M ^Q k ^ e^'(k x e,(k)). 
We set 

(8) 2(a) := Zf(a)+Zf(a), 2(a,6) := U e i(8)Z(a)U e i(0y\ 

Z(8) := 2(0,6), and Z := Z(0,0). 

We consider the Feshbach operator J 7 p^- ) (H g (9) — z) as an operator on Ran P(6>). 
Similarly, we consider 2(a) := Zj(a) + Z° d (a) and 2 (a, 8) as operators on 
RanP e ij and RanP e jj((9) respectively. 

We are now able to formulate our main result. It will be proven in Section H 

Theorem 1. Let < e < 1/3 and g small enough. Let <f>i and 4>2 be normalized 
eigenvectors of H e i with eigenvalue Ej and <£>i := fa ® Q. Assume moreover that 
the imaginary part ImZ := -^{Z — Z*) of Z is strictly positive on RanP e ;j. Then, 
in terms of a dimensionless time parameter s > 0, 

(*i,e-" H «3 2 ) = (fa,e- is ^-^ z ^ 2 ) + b(g,s), 

where \b(g, s)\ < Cg e for some C > 0. 

The theorem has the following immediate Corollary: 

Corollary 2. Under the assumptions of Theorem^ if < r :— g 2 s is kept fixed, 
if 4> := <pi — cf>2 are eigenvectors of Z with eigenvalue T, and if $ := <fi ® CI, then 

lira | {^e -4 "^*)! = e - Tlmr . 
gio 

We close the introductory section with the following remarks: 

Remark 1. The theorem can be rewritten in terms of the original operators: 
Let <p[ and <j>' 2 be normalized eigenvectors of H' el with eigenvalue 2a 2 mc 2 Ej and 
:= # ® fi. Then e^** -1 *^) = e-" 32 ^^-^')^) + C(a 3e / 2 ) = 

(01, e-^^^^-s 2 ^^) + C(a 3e / 2 ), where & <g> = CT[0J <g> fi]. 
Here 

^2 



(9) Z':= 



8a 4 m 3 c 2 



lim / dk K A ^tj-. Pit jeJk) ■ p'P'i ,■ 



x [PeljKl ~ ^ 2 mc 2 E 3 + hc\k\ - iel^P^e^k) -p'P^ 



W dfc *'(|fc|) 2 , , 

- A eR3 fepfe] 4^fej ^J C "W ' P ^U^W ' P feX,i 



P' el i is the projection onto the eigenspace of H' el belonging to the eigenvalue 
2a 2 mc 2 E i: p'j := -ihV X] and p' := YT =iP'y 
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Remark 2. Note that the matrix Z{a) depends on the fine structure constant a, 
since the coupling functions defined in equations ([5|), (J6j, and ([7|) do. Thus, due to 
the exponential decay of the eigenfunctions of the electronic operator, Z(a) can be 
developed in a power series in a = g 2 ^ 3 . The zero order term corresponds to electric 
dipolc (El) transitions, the higher order terms to magnetic dipole transitions as 
well as to higher order electric and magnetic transitions. We have for some C > 

(10) g 2 \\Z{a,9) - Z{9)\\ <Cg 2+2 '\ 

It is easy to see that the imaginary part of Z is (see also O Formula (IV. 19)]) 

j'-i 



(11) lmZ = ir^2Y^ [ duiEi-Ej) 2 

_-_n i o J |uj| = 1 



=0 li=l,2 



Ak(E — E-) 2 

X 47r 2(^. _ ^ djl^M • P]Pel,iK(.v) -P\ P el,j 



r^(£ 3 - £*)k(.Ej - EifP e i tj pP eUi pP ehj: 

6 »=0 



In the last step we used the relationships J2u=i 2( e A*( w ))wi( e A'( to '))« = fim,n — u m u> 



and J du)U! m u n = - 3 "''" , where <5 mj „ is the Kronecker symbol. Moreover, p :— 
^2j—iPj and the expression P e ijpP e i^pP e ij indicates a Euclidean inner product. 
Analogously, we set x := XwLi ^i" Using the commutation relation 

[x, H ei } = 2ip 

we find 

2 

(12) ImZ = 3 23(^3 - E ^) 3 < E 3 - EifPeljxP^iXPelJ. 

i=0 

We analyze equation (fT2")l for the case of a hydrogen atom in Appendix B. We 
show there that ImZ is indeed strictly positive unless j = 1. If j = 1, ImZ has 
a zero eigenvalue, since the 2s state of hydrogen cannot decay via electric dipole 
transitions. However, the 2p states can decay via an electric dipole transition. It 
would be interesting to prove time decay estimates also in the latter case. 

Note that the transition rate is proportional to g 2 a 2 oc a 5 , in accordance 
with physics textbooks (see e.g. [5J Section 59]). 

Remark 3. The eigenvectors of H e i are analytic vectors for the generator of dila- 
tions, and therefore U e i{9) : ker(ff e ;(0) — Ej) — > kcr(H e i(9) — Ej) is a (bounded 
and bounded invertible) mapping between finite dimensional vector spaces. (The 
latter is true for |Im(9| < ir/2.) This implies that the matrices Z(0) and Z(9) are 
similar. In particular, the bounded operators [— g 2 Z(0)<8)l/ + e~ e l e ;<g)i?/]| Ran p(o) 
and [— g 2 Z{9) <£> 1/ + e~ e l e i ® Hf] | Ran p iq\ are similar (cf. [H Section 3]). This fact 
will be used in the proof of Theorem [1] and in Section [3] 
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2. Proof of the Main Result 



In this section we prove our main result. The technical estimates needed in 
the proof are collected in a series of lemmas and deferred to Section [3J For the 
proof we need the operator (see [3l Formula (IV. 67)]) 



(13) QW(*):=V / dkP(0)[w^J(k^)®l f ] 

P(6)(\k\) 



11=1,2 



H el {6)+e-^(H f + \k\ 



defined on RanP(#) and for z G A(8, e). Here we used the definition P(8)(\k\) := 

Pei.j{&) <8> 1/ + P e i,j(0) <8> XH f + \k\>p a - Moreover, we need the the operator Q^\z) 

on RanP^f?), defined by [Q ( o 0) (z)<f>]®tl := Q {9) (z)[<f>®Q] for all <j) G RanP eiJ (6»). 
It is defined by the formula 



(14) Q$\z) 

Y] / dkx\k\>p o Pei,j(9)w { a l(k,fi) 



+ V / dkP elJ (e)w^l(k,fi) 



Ej +e- vd \k\ - z 

Pel AO) 



H el (9) + e-^\k\ - z 



w[ 6) (k,») P eld (9) 

W§(k,ri Pel AO)- 



We remark that both operators are analytic for z G A(S,e). This follows from the 
fact that the resolvents in their definitions can be bounded uniformly in z G A{8, e) ■ 
(See the proof of Lemma[3]for a proof in the case of Q ( o\z), The proof for Q^(z) 
is similar and uses additionally the spectral theorem for Hf.) 

Note that by assumption there exists a constant c > such that Im Z > c. 
Since Z is bounded, there are constants a, b > such that NumRan Z is localized 
as NumRan Z C A(c, a, b), where A(c, a, b) :— ic + [~a, a] + i[0, b] (see Figure [5]). 
We set v := min{$, arctan(c/(2a))}. 

Finally for w G C and r > we define D(w, r) :— {z E C\\z — w\ < r}, and 
for A C C we set D(A,r) := {z G C|dist(z,A) < r}. The notation [z,w] denotes 
either the line segment between z G C and to € C or a linear contour from z G C 
to w € C. Accordingly, [z\, w{\ + [z2, W2] is to be understood either as the sum of 
the sets \z-\_, w\] c C and [Z2, W2] C C or as a generalized contour. 



Proof of Theorem^ First, we show that we can introduce a spectral cutoff with 
an error of O(g): We choose a function F G C§°((Ej - 8/2, Ej + 8/2)) with < 
F(x) < 1 for all x G [E ] -8/2,E j +8/2] andP(a;) = 1 for all x G [£*.,— 5/4, 
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By the almost analytic functional calculus [HI E] 
F(H g ) = - f dxdy^fiiHg-z)- 1 

TT J OZ 

= - f dxdy^P-(H - z)- 1 - - [ dxdy^M(H g - z)- x W g {H Q - z)" 1 

TT J OZ TT J OZ 

where F G (C) is an almost analytic extension of F(x) with | \ = O ( | Im z \ 2 ) . 
Since 

I f dxdy^{e isH <>$i,(Ho-z)- 1 $2) = (e isH °$i,$2) 

TT J OZ 

and 

|- / dxdy(e UH '* u ^(H g - zy'WgiHo - z)-^ 2 )\ 

TT J OZ 

< - [ (fadvii^iiii^nK^ - *)- i ni(^ - zr'iw^w < Cg, 

tt J Oz 

we find that ($i, e - ls ^$ 2 ) = ($1, e~ %sH <> F(H g )® 2 ) + 0(g). 
Analogous to [15| and [4] , we can write 

($ 1 , e -^F(JJ 9 )$ 2 ) = --^Um f d\e- iXs F(\){f(0,\~ie)-f(0,\ + ie)] 

ZTTl elO J 

= -^~J dXe-^F(X)[f(9,X)~f(e,X)], 

where f(9, X) := (^(fl), jjjy^ ^)) with ^(0) 0,(0) <g> fl and &(0) := 
U e i(9)(f>i. We used Stone's theorem in the first step. In the second step we used 
the analyticity of H g (6) and the fact that H g (9) has no spectrum in the interval 
[Ej - 5/2, Ej + 6/2} (see H Theorem 3.2] and also Corollary E below). 

Noting that (^ 1 (9), 7T -^ J - I ^ 2 (9)) = T P(e) {H g {9) - A)-^ 2 (0)) (see 

3, Formula (IV). 14] and also Lemma |A.6|) and using the resolvent equation, we 
obtain 

f(9,X) = (M0),Fp(8)(H g (9) - X^MO)) 

= (MO), [Ej - A - gpQi'H^MO)) ~ (MO), [Ej - A - g'Q^iX)}- 1 <g> 1/ 
x [F P{e) {H g {9) - A) - (Ej - X + e- l el <g> ^ - 3 2 QW(A))P(0)] 

x T P{e) (H g (9) - X)- l M8)) =■ M A ) + A), 

where /(#, A) is the first term in the sum. The strategy is now to move the contour 
for the first term in order to pick up a pole contribution (see Figure [1]) , and to 
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estimate the second term on the real axis: 
tl\e~ iXs F(X) [/(#, A)-/(0,A)] 

d\e~ iXs F(\)[B(9 , A) - B(9, A)] + J dze~ lzs F(z)[f(6, z) - f(9, z)] 

O1+C5 

+ J dze~ lzs [f(9, z )-f(8,z)]- J dze^ zs [f(9, z)-f(9,z)], 

C2 + C3+C4 Co 

where we set C := C a + C 2 + C 3 + C 4 + C 5 , with C x := - <5/2, Ej - 5/4], 
C 2 := [E j -5/4,E j -5/4-ig 2 -*/2],C 3 := -(5/4-z ff 2 -72, ^- + 5/4 - i 5 2 - e /2] , 
d := [Ej + 6 /4- ig 2 - £ /2, E j + 5/4} and C 5 := [E 3 + 8/4, E 3 + 6/2} . C is a suitable 
contour to pick up the pole contribution from f{9,z). The analyticity properties 
required for this process will be discussed below. Note that the contour C cannot 

(9) 

simply be moved down much further, since Q (z) may have singularities outside 
of A(6,e). 



E- 



5/2 E-6/4 



? — H 



E ; +5/4 E ; +5/2 



C, 



"I 



'g 2as 



E,-gA(c,a,b) 



D(E^/2) 



C 5 



c 4 



D(E r g 2 A(c,a,b), Cg 2 >W 



Figure 1 . The integration contour 

Estimates on the real axis: We divide the integration interval [Ej — 6/2, Ej + 
5/2} into two parts: On [Ej — 5/2, Ej +5/2} \ (Ej — pa/2, Ej + po/2) we use Lemma 

I A. 71 and Lemma [7] to obtain |-B(0, A)| < C ■ ( s i l - l $y(\\_E-\-Cg 2 )' 2 • ^ mce 

oo oo 
( ■ <n-2 [ dX I ■ a\-2 -2 f dX 

(sin-i?) / -r = (smi?) g 



, (X~Cg 2 ) 2 v ' » J (A-C7) 2 

2 2 

= (shr^V 2 -^-^ = 0(^ 2 5 - 2 +^), 

we see that the error term for this region is of the order d~ 2 g 3e . On (Ej — po/2, Ej 
po/2) we estimate using Lemma lA.7l and LemmalU \B(6, A)| < C-Q~ 2 - 



2+e 
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Since J dX ^■- 9 xy^+c i g i ^ s eas ily seen to be of order g e , and the same estimates hold 
for B(9, A), the estimate on the real axis is proven. 

Estimates on the contour C: We estimate the integral J c \e~ tsz \\f(9, z) — 
f(6,z)\\dz\. Note that 



1 



- (0! (9) , 2 (9)} + g 2 (0! (9) , Q { e) (z) — -^-6 2 (9)). 



Thus, the zero order terms of f(9, z) and f(9, z) cancel each other, and it suffices to 
show that the higher order terms are at least of order g e . Since Q (z) is uniformly 
bounded in z 6 A(8, e) by Lemma [3j we estimate using Corollary [4] (see Figure [3]) 



1 -M0))\<O' 



E j -(X-ig^)-g 2 Q^ ) (X-ig 2 -^ (Ej - X) 2 + (g 2 ^) 2 ' 

Thus the integral along C3 of the above expression is easily seen to be of order 
g e . The integral over the remaining contour is of order g 2 , since dist(z, Ej) can be 
estimated independently of g along this part of the contour. The integral of f(9, z) 
can be estimated in the same way. 

Estimates on the pole term: Since Q (z) is uniformly bounded for z S A(8, e) 
by LemmaO the function f(9, z) has no poles in A(S, e)\D(Ej, po/2). It follows by 
Lemma |()] that Ej — z — g 2 Q^(z) is bounded invertible if z G D(Ej, p Q /2) \ (Ej — 
g 2 A(c, a, b) + D(0, C x ■ g 2+e )) C D(Ej,p /2) \ [NumRan(£ J - g 2 Z(0)) + D(0, C x ■ 
g 2+e % i.e., all poles of f(9, z) are in the set W := E ] ~g 2 A(c 1 a, b)+D(0, C\ -g 2+e ). 

Moreover, by Lemma [6] we have the estimate \\(Ej — z — g 2 Qo (z))^ 1 ]] < 
Cdist(z, NumRan(£y - g 2 Z))~ 1 for some C > if z e D(Ej,p /2) \ W. In order 
to estimate the pole terms, we choose a contour Co around W such that the length 
of the contour and its distance to W are of order g 2 . A possible choice is Co = 
[E j +g 2 (-(a+c/2)-ic/2),E ] +g 2 (+(a+c/2)-ic/2)] + [E :j +g 2 ((a+c/2)-ic/2),E j + 
g 2 {+(a + c/2) - i(b + 3c/2))] + [Ej + g 2 {+(a + c/2) - i(b + 3c/2)), Ej + g 2 (-{a + 
c/2)~i(b + 3c/2))} + [E J +g 2 (-(a + c/2)-i(b + 3c/2)),E J +g 2 (-(a + c/2)-i(c/2))}. 
We now use the expansion 

(MO), (Ej - z - g 2 Q^(z))- 1 M0)) = (MO), (Ej - z - g 2 Z(9))- 1 M0)) 
+ g 2 (M0), (Ej z g 2 Q i 6) (z)r 1 (Q { 6) (z) - Z(9))(Ej - z - g 2 Z (9))- 1 M®)) ■ 



The integral over the first term gives the claimed leading term, the second term is 
of order g e by Corollary I A . 9 1 and Lemma [51 
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Since by the above considerations (with 8 replaced by 9) the function f(9, z) 
has no poles in the lower half-plane, there is no pole contribution from this function 
(see also Remark 2]). □ 



3. Technical Lemmas 

As in [4], we need estimates on the numerical range and on the norm of 
the inverse of various operators. We make use of numerous results shown by Bach, 
Frohlich, and Sigal 4J , which are summarized in Appendix A. We use the following 
definitions from [4]: For r\ > such that Ej + 5/2 < £ — n we define Pdisc(0) := 

Ei-.EtKX-r, P i( ) and Pdisc(0) ■= 1 - P«K.c(0)- 

Since the operator valued function is relevant for the location of the pole 
term in the time decay estimates, we need certain properties: 

Lemma 3. Let d sufficiently small and g e / sind < 1/2. Then Qq (z) is uniformly 
bounded for z £ A(S, e) . 

Proof. The proof follows [3l Chapter IV], using, however, the following estimates: 
For the first summand in equation (|14ll . we use the estimate |e _l,? |fc| + Ej — z\ > 
|Im(e _ ^|jfc|+£y -z)\ > |sini?|fc| - g 2 ~ £ \ > | sin. i?| • ||fc| - p /2\ > l/2sini?|fc|, for 
\k\ > po- For the second summand in (flT]) , observe that for all Ei with i ^ j we 
have \Ei + e' e \k\ - z\ > sini9(5/2 - g 2 ~ e > 1/45 sin and that by Lemma [A~T1 



(H el (0) - (z - e-™\k\))- L P disc (0) < 



S — r\ — Re z + cos ^|fc| 

□ 



This has the following immediate corollary: 

Corollary 4. There exists a constant C > such that for all z G A(S, e) 

NumRan(^ - g 2 Q^\z)) C D(Ej,C ■ g 2 ). 

We use the following lemma to estimate the inverse of the Feshbach operator: 

Lemma 5. Suppose A is a bounded operator on a Banach space and let A' be 
similar to A, i.e., there exists a bounded, bounded invertible operator G such that 
A' = GAG^ 1 . Moreover, let B be a another bounded operator. Then for any q > 1 
and for all z ^ £)(NumRan(yl), q ■ ||£>||||G|| • ||G _1 ||) the following estimate holds: 

\\(A' + B - z)-^ < \\G\\ ■ IIG" 1 !!— — ■ dist^NumRanM))- 1 . 

q-l 

In particular, a(A' + B) C D(NumRan(A), q ■ \\B\\\\G\\ ■ || C~ 1 1| ) . 

Proof. First, observe that for all z (/ NumRan(A) by similarity 

\\{A' ~z)- 1 \\ < \\G\\ ■ KG" 1 !! ■ IKA-z)- 1 !! < ||G|| • HG- 1 !! •dist(z,NumRan(A))- 1 . 



12 



D. HASLER, I. HERBST, AND M. HUBER 



By a series expansion we obtain for z £ D(NumRan(^4), q ■ \\B\\ ■ ||G|| • ||G 1 ||) 

oo 

{A' + B-Z)- 1 = (A' - z)- 1 ^[-B(A l - z)- 1 }" . 

n=0 

Taking the norm of both sides implies the claim. □ 

Following [4], we control the Feshbach operator !F P ^(Hg(9) — z) for z G 
D(E 3 ,p a /2) as follows (see Figure [2]): 

Lemma 6. Let < $ < 0q and < g <C 1? small enough. Then the following 
statements hold: 

a) There are constants C\,C% > such that J r P rg\(H g (6) — z) is bounded 
invertible for all z G D(E 3 ,p /2)\D(NumRan(Ej-g 2 Z(0)(g)lf + e~ e l e i(g) 
-ff/)lRanP(o), Ci • g 2+e ), and for A G [Ej — po/2, E 3 + p /2] the estimate 

\\T P{e) (H g (6) - A)- 1 )) < C2 = 

smi/y {Ej — \y + cg^ 

holds. The same holds for (Ej — z — g 2 Q^ B \z) ® If + e~ s l e i )|R an p(e) ■ 

b) There is a constant C > such that for all z G C \ NumRan(£j — 
5 2 ^'(0))|RanP c! , i (o) th e operator (Ej - z - g 2 Z(6))\ KanPcl ](g) is bounded 
invertible and fulfills the estimate 

(15) \\[(E ] -z-g 2 Z(e))\ RanPclj{9) }- 1 \\ 

c 



< 



dist(£,NumRan(.Ej - g 2 Z(0))\ RanPcl 

There are constants C\,C2 > such that for all z G D(Ej, po/2) \ 
D(NumRan(i?j — g 2 Z (0))\ Ran Pcl ^ ), C\ ■ g 2+t ) the operator (Ej — z — 

9 2 Q < o\ z ))\'RanP el j(e) is bounded invertible and fulfills 



(16) \\[(E J -z-g 2 Q^(z))\ K&aPsli(e) r l \\ 

° 2 



dist^NumRan^ - g 2 Z(0))\ RanPal j{0) )' 

Proof. By similarity (cf. Remark [3]) we obtain immediately for some C3 > that 

|| [(Ej - z - g 2 Z(6) ®l f + e~ e l el ® H f )\ Kan p^'W 

< Ci •dist(z,NumRan( J B i - g 2 Z(0) <g> 1/ + e~ e l e ; <g) H f )\ Ran p^y 1 . 

By Lemma [A. 71 Corollarv lA.91 and Lemma [5] there are constants C\, C2 > such 
that 

(17) \\F P{e) (H a (0) ~z)- l \\ 

< C 2 dist(z, NumRan(Ej - g 2 Z(0) 1/ + e~ e l e i ® i?/)|R a nP(o)) _1 
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E-5/4 



/ 



E r g 2 a E^a 



£,+5/4 




i X t i 

I eg 



Figure 2. The numerical ranges of the operators Ej — g 2 Z(0) 
and NumRan^ - g 2 Z(0) <g> 1/ + e- Q \ el ® #/)| Ran p(o) 

follows for z ^ Z?(NumRan(^- - g 2 Z(0) ® 1/ + e - e l e; <g> ff/)| Ran p(o)), d • c/ 2+e ). 
It follows that NumRan[(-6f 2 Z(0)«il/ + e- e l ei (g)iJ / )| Ran p( )] C ~5 2 A(c, a, 6) + 
e~ [0,po] (see Figure [5]). By geometrical considerations, we see that this set is 
contained in the conical region — i| g 2 — i{re l ^\ — (v — \) < 4> <v — |,r € [0, oo)}. 
This, in turn, implies the claim. The claims in b) follow by the same reasoning. □ 

We do not see that the estimate of Lemma [6] a) is true for A S [Ej —6/2,Ej + 
6/2] \ (Ej - p /2,Ej + pa/2) as used in [H Proof of Theorem 3.5] (see also the 
remark after Lemma [A.8|) . Thus we bound !Fp(a\(H g (6) — A) -1 differently in that 
region in the next lemma. 

Lemma 7. Let < d < 9q and < g <C $ small enough. Then J 7 p^(H g (9) — z) 
is bounded invertible for all z S A(6, e) \D(Ej, pa/2) and there is a constant C > 
such that for g small enough and z € A(5, e) \ D(Ej, po/2) the numerical range of 
Tp(Q}{Hg(ff) — z) is localized as 



NumR&n(F P{e) (H g (9) -z)+z)C D(E, + e~ e [Q, p ], Cg 2 ). 
In particular, for A 6 [Ej — 6/2, Ej + 6/2} \ (Ej — po/2, Ej + po/2) the estimate 



Proof of Lemma^ We have by Lemma [O that \\P(0)Wg(9)P(0)\\ = 0(g 2+t ). 
Therefore, it suffices to show that 




\\p(e)w g (9)p(0)[P(e)(H g (e) - z)p(e)r 1 p(e)w g (e)p(e)\\ = o( 9 2 ). 
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E-§/2 E r 8/4 



C, 



C, 




£,+5/4 E.+8/2 



D(E j V[0,p ] > Cg 2 )- 



Figure 3. Global localization of the numerical range of Ej 
g 2 Qo(z) and of the Feshbach operator f r p^(H g (9) — z) + z. 



Following 4, Proof of Lemma 3.14], we use a Neumann expansion: 
(18) P(9){P(9)(H g (6) - z^)}- 1 ^) 



3C 



= ]T P(9)[P(9)(H (9) - z)P(9)]-ip(6) 



71=0 

X 



\-P(9)W g (9)P(9) [P(9)(H (9) - z)P{9)] 1 P{9) 



This expansion is valid if z £ A(S, e) with Im z > C for some C > (independent 



of g.) We define B g (p) := H e i{9) ® 1/ - Ej 



(l et ® H f +p) as in @]. The right 



n=0 



handside of equation ([18]) is equal to 

oc 

(19) ^ ii? e (po)r i/2 ii? e (po)r /2 ^) [^(f?)(^ (^) - ^)^)]- i ^ ) | i? _ (/0o) |i/2 

-|S e -(po)|- 1/2 W ff (^)| J B e (po)|- 1/2 |S e (po)| 1/2 PW 

x [P(^)(ifoW-^)PW] PW|S e -(po)| 1/2 ] |S e -(p )|- 1/2 =:^). 

for all z £ A(S, e) with Im z > C. By Lemma IA.4I and Corollary IA.3I the series 
in equation (|19p converges uniformly for z £ A{5, e) and is thus a holomorphic 
function of z £ A(S, e) . Thus 

P(9)[P(9)(H g (9) - z)P(9)]^P{9) = R{z) 

for all z £ A(S, e) by holomorphic continuation and for all z £ A(S, e) 

(20) P(9)W g (9)P(9) [P(9)(H g (9) - z)P(d)]~ 1 P(d)W g (9)P{9) 

= P{9)\B s { Po )\ 1 / 2 \B s {p G )\- 1 l*W g {9)R{z^ 

Note that \\\B e ( Po )\P(9)\\ - \\B e (p )P(9)\\ and \\P(9)\B § ( Po )\\\ = \\P(9)B e (p )l 
Thus, using ||P(0)|i? ff (A))| 1/3 || < \\P(9)\B s (p )\\\ ■ \\\B g (p )\-^\\ - 0(pl /2 ), and 
counting the powers of po in (12"01) , the first claim follows. The estimate on the 
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inverse follows by geometrical considerations. The claim on Ej — z — g 2 Q^\z) 



Note that due to the appearance of the interaction W g (9) on both sides of the 
resolvent [P(9)(H g (9) — z)P(9)]~ 1 and due to the projections P{9), the divergence 
of the resolvent for po — > is completely eliminated (see also the remark after 
Lemma IA.6|) . 

We use the following corollary instead of [H Theorem 3.2]. 
Corollary 8. Let < d < 9q and < g <C small enough. Then 



for some C > 0. In particular, the interval [Ej — 8/2, Ej + 8/2] is contained in the 
resolvent set p(H g (9)). 

Proof. By Lemma J 7 p/g\(H g (8) — z) is bounded invertible for all z S A(8, e) \ 
D(Ej,p /2). By Lemma [51 it is bounded invertible for all z E D(Ej, p /2) \ 
D(NumR&n(Ej-g 2 Z(0)(g)lf + e~ e l e i(g)Hf)\ Ran p( ),Ci-g 2+e ). Lemma [Ql implies 
the claim. □ 

Appendix A. Estimates Taken from Bach, Frohlich, and Sigal [3] 

In the appendix, we quote some important technical lemmas from [1] , which 
we frequently use. We do not give their proofs, since they are very lengthy. However, 
for the orientation of the reader, we describe the essential points of the proofs in 
words. 

A.l. Existence of the Feshbach Operator. First we need certain relative 
bounds on the interaction and bounds on the resolvent. 

Lemma A.l (4 ], Lemma 3.8). Let z G C with Re z < E - 77. Then, for |0|(1 + 
(S — rj — Rcz) -1 ) sufficiently small, H e i(6) — z is invertible on Ran Pdisc{0) and 



This lemma is proved by using that the estimate holds for 9 = with constant one 
(instead of two) and using that H e i(0) — H e i(6) is relatively H e i(Q) bounded. We 
remind the reader that as in |4 we define Bg(p) := H e i(9) ® If — Ej + e~ e (l e ; CS> 
Hf + p). Note that A(8, e) c p(P{9)H (9)). 

Lemma A. 2 ([4], Lemma 3.11). There exists a constant C > such that for 

— 1/2 

< ■& < 9q, for all g with < gp Q < 1/3 and < po < (8/3) siwd, and for all 
z e A{8, e) 



follows from Lemma [31 



□ 



A(8, e) \ (Ej - D(g 2 A(c, a,b),C ■ g 2+e ) + e^O, p ]) C p(H g (9)) 



(P dlsc (9)H el (9)P dlsc (8) - z)- x P disc (6)\\ < 2(S - r? - Re z) 



(21) 



\\Be( P0 ) 



P(8) 




H (9) 



— z 
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The proof of Lemma IA.2I is based on Lemma IA.1| the fact that H e i (9) restricted 
to Pdisc{9) is similar to a self-adjoint operator, and various other estimates on the 
resolvent of H e i(9) as well as the application of the spectral theorem for Hf. 
The following Corollary was used in [2] : 

Corollary A. 3. There exists a constant C > such that for < $ < 9 , all g 

with < gp < 1/3 and < po < (5/3) sini?, and for all z G 4(5, e) 



I^Mr /2 Tr!^l%(po)l 1/2 ll 



c 

H (9)-z^^n n ^j- 



Proof. By taking adjoints in equation (j2"Tj) . we find that II Ho(6)-z \B§(Po)\ II — %■ 
The claim follows by complex interpolation. □ 

Lemma A. 4 ( 4J, Lemma 3.13). There is a constant C > such that for < $ < 
9q sufficiently small, #i, #2 € {±i$} and for all p > 



c , 



(22) iii^^i- 1 / 2 ^^)!^^)!- 1 / 2 !! < 



p 



-1/2% 



< 



The proof of Lemma [O] uses that ||4 K (x)_Vll < C||#/ V>ll and ||4 K (a;) + ?A| 
C||(if/ + 1) 1/2 ?A|| for some < C and all ^ in the domain of H 1 / 2 , that H e i(0) - 
H e i(&) is relatively H e i(0) bounded, and some other estimates. The term propor- 
tional to p~ x l 2 is due to the +1 in the bound for the creation operator A K (x) + 
and to the appearance of a similar constant in the estimate for the relative bound- 
edness of H e i(0) — H e i(9). Note that the symmetric form of the estimate (|2"2")l is 
essential. Estimates on W g (9)\Bg 2 (p)| _1 lead to worse behavior as p — > 0. 

Lemma A. 5 ([4., Lemma 3.14). There is a C > such that for § € (0, 9o), po < 
(5/3) sintf , < .gpo 1/2 < i? 2 ; and /or a// z G .4(5, e) the operator P(9)H g P(9) - z 
is invertible on RanP(#) and fulfills 

\\^(e)H g (6)P(d)-z]- l P(9)\\<^-. 

The proof of Lemma IA.5I uses Corollary IA.31 Lemma IA.4I and a Neumann series 
expansion. 

Lemma A. 6 ( 4 , Lemma 3.15). Assume that $ G (0, 9q). Let po < (5/3) sm-d and 
< gp <C i? 2 . Then for all z G 4(5, e) the Feshbach operator J 7 p^(H g (9) — z) 
defined in equation ^ exists. If z G 4(5, e), then H g (9) — z is bounded invertible 
if and only if the Feshbach operator J- P ^(H g (9) — z) is bounded invertible, and 
the equation 

(23) (H g (9) - z)- 1 = [P{9) - P{9){P{9)H g {9)P{9) - z)- 1 P(9)W g (9)P(9)} 
x T P{e) (H g {9) - z)- 1 [P(0) - P(9)W g (9)P(9)(P(9)H g (9)P(9) - z)- 1 ^)} 

+ P(9)[P(9)H g (9)P(9) - z}- 1 ^) 
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holds, where the left side exists if and only if the right side exists. Moreover, there 
is a constant C > 0, independent of g and 9, such that for all z 6 A(S, e) 

(24) \\{P{9)H g {9)P{9)-z)- l P{9)W g {6)P{9)\\< ^ 



Q 1/2 



and 

Cg 



(25) \\P(9)W g (9)P(0)(P(9)H g (9)P(9) - z)-^ < 



a 1/2 ' 



Equations (|24|) and (|25p are proved similarly as Lemma [A.51 Together with Lemma 
IA.5l thev imply the existence of the Feshbach operator and the validity of Equation 
(|23| (see Theorem IV. 1]). Note that the operator W g {9) in Formulas (JMJ) and 
(f25|) reduces the divergence as po — > in comparison to Lemma IA.5I 



A. 2. Approximations of the Feshbach Operator. The following lemma gives 
an approximation of the Feshbach operator globally for all z G A(5, e) (see [1], 
Lemma 3.16, estimates on Remo through Rem^): 

Lemma A. 7. Let < e < 1/3 and < d < 9q. Then there is a constant C > 
such that for all g > sufficiently small with po < (<5/3) sind, and for all z S 

|| F> (fl )(# fl (0) - z) - (Ej - z + e- l el ®H f - g 2 Q^(z))]P(9)\\ < Cg 2+i . 
Moreover \\P(9)W g {9)P(9)\\ = 0(g 2+e ). 

The lengthy and technical proof of Lemma I A. 71 is based on a Neumann series 
expansion, estimates similar to Lemma [A.41 and the pull-through formula. 

For z sufficiently close to Ej, Q^'(z) can be approximated by Z(a,0) (see 
[H Lemma 3.16, Estimates on Rem.4 and Rem^\). 

Lemma A. 8. Let < e < 1/3 and < d < 9q. Then there is a constant C > 
such that for all g > sufficiently small with po < ((5/3) sind, and for all z € 
D(E 6 , po/2) 

5 2 ||QW(z)-Z(a,0)|| <C 5 2+e . 

The proof requires some additional estimates to eliminate the z-dependence 
of QW(z). However, we not see that Lemma IA.8I holds for all z € -4(<5, e) , which 
seems to be used in |J. 

Lemma [A. 81 and Equation (| 10[) imply 

Corollary A. 9. Under the assumptions of Lemma \A.8\ 

g 2 \\Q^(z)-Z(9)\\<Cg 2+i . 

Remark 4. Note that in order to approximate the Feshbach operator J- p ^(H g (9) — 
z) for 9 — id with i? > 0, the — ie in definition ([3|) has to be replaced by +ie. In 
particular, when considering the spectral analysis of this operator, the localization 
of the numerical range and of the spectrum have to be reflected about the real 
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Appendix B. The Hydrogen Atom 

In this section we discuss the applicability of the presented method to the 
hydrogen atom. In particular, we show that Im Z is strictly positive unless j = 
1. For compatibility with physics literature, we number the eigenvalues of the 
hydrogen atom according to the principal quantum number n — i + 1 . We denote 
the corresponding eigenvalues by (£„, i.e., <£„ = Ei for all i > 0. We will ignore the 
(trivial) spin dependence of Z = Z(0, 0) in this appendix. 

B.l. The Hydrogen Eigenfunctions. We define the associated Laguerre poly- 
nomials (see [6J Formula (3.5)]) for A, /j 6 No with < [i < A by 

and set (see [SI Formula (3.16)]) 

< 26 > : = -jg ^X^ "'-"™ (;)' WW")- 

Note that the Hamiltonian in [B] has an additional factor of 1/2 in front of the 
Laplacian, so that the radial functions and certain other quantities have to be 
adapted accordingly. We would like to warn the reader that there are different 
conventions for the indices of the associated Laguerre functions. 

For n £ N and I, m £ Z with < / < n ~ 1 and —I < m < I the normalized 
eigenfunctions to the eigenvalue <£ n are 

(27) u„,i, m (r, 9, 0) := Rnj(r)Y ltm (p, <p), 

where the Y^ m are spherical harmonics (see [6l Section 1]) and we introduced polar 
coordinates by 

x = r sin 9 cos <j> 
y = r sin 9 sin <j) 
z = r cos 9 

with < 9 < 7r and < <f> < 2ir. Note that in this appendix x, y, and z denote 
the cartesian coordinates of the electron, contrary to the main part of the paper, 
where x and z have different meanings. Moreover, note that the eigenvalues (£„ 
are n 2 -fold degenerate. 

B.2. Selection Rules for Dipole Transitions. In this subsection, we give some 
important results from [6\. We define 

oo 

(28) R^f := J drr 3 R n ,^(r)R nA (r). 
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These integrals have been evaluated by Gordon [13] (see also Section 63]). 
Below, we need (see (H Formula (63.4)]) 



(29) |i?2,i°l = 2- 



'2 15 n 9 (n- 2) 2n 
3(n + 2) 2 ™+ 6 



For the dipole moments (it n ',i' m'j z u n ,i,m) one finds (see [H Formula (60.11)]) for 
all n, n' S N that 

(30) (v,l',m', zun^i^m) = unless 2 = 2 ± 1 and to' = m. 
Moreover, we will need the relation 

(31) Ki',o,o, zu 2 ,i,a) = J^ R 2,i°- 

The selection rules given in [SJ Formula (60.11)] imply immediately 

(32) (Un' .1' ,m' i XUnJ,m) — (^n' ,V ,m' ; yUn,l,m) — 

unless I' - - 1 ± 1 and m' = m ± 1 

B.3. The Imaginary Part of Z. In this subsection, we show that the method 
presented in this paper applies to the hydrogen atom, except for the case n = 2. 

Theorem B.l. Fix n € N and consider 
1 



X 



[P e ljxP e l^iXP e lj + P e l j y P e l P e l J + P e lJ zP e l ti zP e l t j] 



for j = n — 1 as in equation (| 1 2|) . TTien /or a22 2, m, 2', m' S No wi£/i < 2 < n — 1, 
-2 < m < I, < 2' < n - 1, and -2' <m' <V 

unless I = 2' and m = m! , and for all 2, m G No with < I < n — 1, — 2 < m < 2 

(u nt i im ,1m Zun t i, m ) > 
unless n — - 2. 7n particular, ImZ is positive, unless n = 2. 

Proof. Off-diagonal matrix elements: Since ImZ is invariant under rotations, it is 
diagonal in the basis {u rai z jm | < I < n— 1, — / < to < I}. This can also be verified 
using the explicit formulas for the dipole matrix elements in [6l Section 63] . Note 
that the matrices P e i t jxP e i ti xP e ij, P e i^yP e i ti yP e ij, and P e i,jzP e i,iZP e ij are not 
diagonal separately. We would like to mention that also the real part is diagonal 
in the basis {u n ^ tm \ < I < n — 1, — I < m < I}. 

Diagonal matrix elements: Let us first remark that the matrix clement 

("2,0,0, [Pel.lxPelflXPel,! + Pel ,W Pel ,0V Pel .1 + Pel .1 Z P e l ,0 Z P e l ,l]w2 ,0 ,o) 

vanishes by the selection rules (f3"2")l and (|3"0|). 
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Suppose now that n > 3. We have to prove that there is an i < j = n 
such that for all <fi € RanP e jj 



Since ImZ is diagonal in the basis {u„ tm j\l — . . . n— 1, m = —I, . . . , I}, it suffices 
to show 

^ \\Pel,iPvU n ,Lm\\ 2 > 
v=x,y,z 

for all < I < n — 1 and —l<m<l. For the case Z = 0, to = it follows 
from equations (|3Tj) and (|29f that the transition (n, 0, 0) — > (2, 1, 0) is an allowed 
electric dipole transition, since z% \ ' > 0. Consequently (un^o^^o,!)) > 0. 

Thus, it suffices to consider the case / > 0. The proof is by contradiction. 
Assume that Y^ v =x, y ,z \\Pei,iPvUn,i,m\\ 2 = for all i < j = n - 1 and some l,m. 
This would imply that for v — x,y, z 

For I > 0, it is easy to see by Equation (f2"5)) that p v u n ,i,m & Dom(H e i) and, using 
partial integration and the fact that u nt i^ rn (0) = 0, we see that 

^ {,Pv^n,l.m^ H e ip v 1l n i 7n ^ — Ej ^ {jPv^n^l.ra^Pv^n^l^rn)- 
v—x : y 7 z v=x,y,z 

Thus, we conclude by the variational principle that 

HelPv^n,l.m — PjPv^n^ljiri' 

However, 

EjP v U n i rn — HelPv^"n.l,m — ^jPv^7i,l,m ~t~ [-^edPt;]^n,?,m 

for v — x,y,z and 

X 

[H e l,p x ] = -i-j 



with r = \J x 2 + y 2 + z 2 , so that we arrive at a contradiction. 



□ 



B.4. Numerical Illustration. In this subsection we give explicit numerical val- 
ues for the matrix Im Z for the case n = 3 setting the cutoff function k identically 
equal to one. Using Maple and the explicit form of the eigenfunctions in equation 
(|2T|) . we calculate the matrices P e i,oxP e i^ and P e i,\xP e i,2 as well as the corre- 
sponding matrices for the coordinates y and z, where P e ifl is the projection onto 
the groundstate, P e i t x the projection onto the eigenspace belonging to £2, and 
P e i,2 the projection onto the eigenspace belonging to (£3. With these matrices, we 
calculate Im Z according to equation (|12p . The numerical values for other principal 
quantum numbers could be calculated in the same way. 

The matrix ImZ (and also Z) is diagonal in the basis {u3./, m |0 < I < 
2, — I < to < I}. The diagonal elements depend only on I, but not on to. We find 
(u 3 ,o,o,ImZu 3i o,o) = 19 53i 25 , {u 3A . m ,lmZu 3A . m ) = 25 7 n 3 8 4 2 3 00 for -1 < to < 1, 
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and (113,2,771, Ini Zu^^.m) — ^2^125 ^ or — ^ < m <2. Let us remark that the eigen- 
values of 2 • (2a 5 mc 2 /S)Im Z are precisely the inverse lifetimes t~v of the corre- 
sponding eigenstates of the hydrogen atom. The additional factor two is due to the 
fact that lifetimes are defined via survival probabilities and not via survival ampli- 
tudes. Inserting a = 7.29735-10" 3 , m = 9.10939-10~ 31 fc 5 , c = 2.99792-10 8 m/s and 
h = 1.05457 ■ 1CT 34 Js we find t 3 , ,o = 1.58303 ■ 10~ 7 s, r 3 ,i, m = 5.26860 ■ 10~ 9 s for 
— 1 < m < 1, and T^^.m — 1.54593 T0~ 8 s for —2 < m < 2. Experimental values for 
these lifetimes are not very precise. We quote a value of T 3 ,i, m = (5.5±0.2) x 10 _9 s 
given in [S]. H find a value of r 3 ,i, m = (5.58 ± 0.13) x 10" 9 s and PJ find 
T 3, i,m = (5.41 ± 0.18) x 10 _9 s. Notice that the experimental values are in rea- 
sonable agreement with the calculated value. 
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